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Abstract 

We study the evolution with temperature of quasiparticle bound states 
around non-magnetic impurities in dx2_y2-wave superconductors. The associ- 
ated local density of states has a fourfold symmetry which has recently been 
observed in Zn-doped Bi2212 using scanning tunneling microscopy (STM). 
From the corresponding Bogoliubov-de Gennes equation we find that with in- 
creasing temperature the magnitude of the bound state energy increases and 
the amplitude of the fourfold contribution to the spinor wave functions de- 
creases. In the pseudogap regime above Tc the fourfold angular dependence 
of the local tunneling conductance persists as long as the superconducting 
fluctuations are sufficiently strong to support a finite local order parameter. 
Once the gap function vanishes completely, the angular structure of the bound 
state wave function becomes featureless. These effects should be observable 
in STM studies of impurity doped high-temperature superconductors. 
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Introduction: It is well known that impurity doping, such as Zn or Ni substitution in the 
high-temperature superconductors LSCO, YBCO, and BSCCO, is a useful tool in demon- 
strating the underlying nodal structure of their order parameter. It also provides a semi- 
quantitative test of the BCS theory for dx2_y2-wave superconductors. Jlj-Q In particular, an 
examination of the thermodynamic and transport properties in these compounds suggests 
that Zn impurities can be modeled with a scattering potential in the unitary limit. 
On the other hand, little is known experimentally about the local structure of the impurity 
bound states around the Zn sites despite numerous theoretical studies on this question. [^,|^ 
Recently, Pan et al. have provided the first scanning tunneling microscope (STM) images 
of the local tunneling conductance around Zn impurities in Bi2212 at low temperatures, T 
4.3 K, and fixed at iE', where E ^ A/30 is the binding energy of the impurity bound 
state. The corresponding wave function was shown to exhibit a fourfold angular symme- 
try associated with the underlying dx2_y2-wave order parameter. We have found that this 
quasiparticle bound state wave function around impurity sites can be described within the 
formalism of the Bogoliubov - de Gennes (BdG) equations for dx2_y2-wave superconductors. 
Within this framework, the main features of these observations, in particular the char- 
acteristic angular patterns seen in the STM images of Zn-doped BSCCO, were reproduced. 

While this previous work was based on an analysis of the BdG equations at zero temper- 
ature, the thermal evolution of the properties of the impurity bound state is of particular 
interest. Here we will study how the angular STM patterns around impurities may pro- 
vide valuable information regarding the nature of the pseudogap in the temperature regime 
above Tc. We will assume that the origin of the pseudogap is dominated by standard su- 
perconducting fluctuations. |8HlT| This assumption is motivated by recent angle-resolved 
photoemission spectroscopy (ARPES) measurements on underdoped BSCCO which have 
clearly indicated an angular dependence of the pseudogap feature proportional to cos^ (20), 
analogous to |A(k)p in the dx2_y2-wave superconducting phase of this compound. [|I2| In 
the following, we will address the consequences of the pseudogap phenomenon on the bound 
state wave function in a semi-phenomenological manner. In particular, we will investigate 
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how the characteristic fourfold symmetry pattern in the wave function gradually disappears 
as the temperature is increased. 

Bogoliubov - de Gennes equations: Let us start our analysis by examining the BdG 
equations for the spinor wave functions in a dx2-_y2-wave superconductor, given by 

Euir) = (^-^ - ^ - Viv)^ uiv) + ^Aidl - d^Mv), (1) 

Evir) = - (^-^ - /i - V(r)^ ^;(r) + ^A{dl - 9j)«(r), (2) 

where /i is the chemical potential, is the Fermi momentum, and V{r) = Vo5^(r) > is 
an isotropic impurity scattering potential, centered at r = 0. In previous work it was shown 
that the spinor wave functions u{r) and f (r) can be expanded in terms of Bessel functions 
of the first kind, leading to the Ansatz P 

n(r) = A exp (-7r) [Mppr) + v^pj^ippr) cos (40)) , (3) 
v{r) = V2Aaexp {—^r)J2{pFf) cos (20), (4) 

where Ji{pF'r) are Bessel functions of the first kind, a, /3, and 7 are variational parameters, 
and A is a global normalization factor. In this expansion, only the leading-order angle- 
dependent terms have been retained. 

Since the spinor wave function v{y) does not contain an s-wave component, it can be 
eliminated from Eq. (1), yielding 

where the contributions of the kinetic energy and the impurity scattering potential are given 

by 



drr {drexp {—'yr)Ji{pFr)) + {I exp {—'yr)Ji{pFr)/r 
K = 



2m drr {exp {—'yr)Ji{ppr)Y 2m 

Note that the contribution of the impurity potential V does not appear in the denominator 
of the last term in Eq. (5) because the spinor wave function f (r) does not have an s- 
wave component. Furthermore, the kinetic energy contribution reduces to 7^/2m for all 



Ji{pFf) with ^ <^ J9f/7- The two parameters of interest which will be determined in 
the following are the temperature-dependent binding energy E{T) and the coefficient of the 
fourfold symmetric term (3{T). Because u{y) contains two orthogonal angular components, 
and the kinetic energy K is independent of angular momentum, it follows that 

Solving these two equations for E and /?, one obtains 

A2 

13 = , (10) 

AVS = [e^ -K^ + V{K + E)- AV2) [e^ -K^- AV2) . (11) 

The quartic equation for E contains one root which will determine the bound state energy. 

To make further progress, we will assume that (i) the impurity scattering is in the unitary 
limit, Vo ^ Aq, leading to a small binding energy E especially at low temperatures, (ii) the 
energy contributions K and V are independent of temperature in the regime of interest, 
T G [0, 3Tc], and (iii) the local amplitude of the superconducting gap function A(T) is small 
but finite in the pseudogap regime above T^.. 

Superconducting regime: Let us first consider the zero-temperature limit, assuming that 
the impurity scattering potential is in the unitary scattering limit. In this case, — > 0, 
leading to 
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V 



K{K' + A^) + At/8K . (12) 



+ Al/2 

Following the assumption of temperature-independent energy contributions K and V, we 
can set K = oAq, where Aq is the gap amplitude at T = , and a is a proportionality 
constant whose physical range can be inferred from experiments. It then follows that 

l + 8a2(l + a^ 
^ " 4a(l + 2a2) ^^^^ 



4 



Furthermore, it is known that at zero temperature (3 < 0, leading to the characteristic 
fourfold symmetric patterns in the local density of states around impurity sites that have 
been observed in recent STM measurements. It is therefore natural to postulate that 
P < for all temperatures, imposing the constraint: a > 2^'^^^ = 0.5946. Therefore, we will 
only consider the physical parameter regime a > 2~^/^. For example, we can set a = 2/3, 
yielding P = — 9/(17-\/2) = —0.3743 at T=0. We have verified that the finite-temperature 
properties of the bound state that will be discussed in the following do not depend strongly 
on the choice of a. 

Normal state: In the opposite limit, T ^ Tc, the gap function vanishes A^(T — ^ oo) —>■ 0. 
Assuming that K and V are basically unaffected by temperature, this gives 

EiT^o.)^K-V^-^L^^^,. (14) 

Hence, the magnitude of the impurity bound state energy E increases with increasing tem- 
perature. The parameter (3 vanishes in this limit, indicating that the bound state wave 
function loses its fourfold angular symmetry pattern at high temperatures when A^(T) 0. 

Pseudogap regime: In the pseudogap regime the amplitude of the local order parameter, 
|A(T)|, is assumed to be small but finite, and by perturbing about the limit T ^ oo one 
finds to leading order that 

A2(T) (l + 4a2) ^ 2a(l + 2a2)A2(T) , , 

and 

' 2x/2(64a6(l + a2)-(l + 8a2))A2 ^ ' 

Hence the magnitude of the binding energy E{T) decreases as T — * due to the progressive 
opening of the energy gap. 

Provided that the pseudogap above Tc arises mainly due to superconducting fluctuations 



T|, it is straightforward to incorporate its physical consequences by assuming a small 



but finite amplitude of the local order parameter above T^. This can be modeled by 



A' 



(T) = - [aUi - {T/T^ff + v/A^(l-(r/T,)3)4 + c , 



(17) 



where C = 0.027572. The coefficient C is obtained from 



C 



2ttT 



(18) 



4A4((T/r,)3 - 1)2 



ma In (T/T,) ' 



where the right hand side is the spatial average of the fluctuation order parameter, and = 
73(3)f ^/2(27rT)^. The approximate ratio T^/Ep ~ 0.03 in Bi2212 has been deduced from 
the low-temperature behavior of the thermal conductivity. ||13| This interpolation formula 
for A^(T) has already proven successful in the analysis of the excess Dingle temperatures in 
the vortex state of the k — {ET)2 salts. 0,|l^] 

In Fig. 1, the temperature-dependent gap amplitude is plotted along with the energy of 
the bound state and the coefficient of the fourfold symmetry term in the spinor wave func- 
tion u{r). The gap amplitude at zero-temperature has been set equal to unity, Aq = 1. 
At very small temperatures, we observe that A(T — * 0) — > 1, E{T — * 0) — > 0, and 
P{T ^ 0) ^ — (v^(a^ + 1))~^, as expected from the discussion of the T — > limit. In 
the limit T ^ cx), we find that A(T ^ 0) ^ 0, E{T ^ -(1 + 40^) (4a + 8a^)-\ and 
(3{T 0) —>■ 0, implying that the magnitude of the binding energy increases with tem- 
perature, and that the fourfold angular features in the spinor wave function disappear at 
high temperatures. The intermediate-temperature regime of E{T) and /5(T) in Fig. 1 has 
been determined by a numerical solution of the coupled equations (10) and (11). In the 
fiuctuation regime around Tc, the coefficient f3{T Tc) is found to be appreciable. There- 
fore, remnants of the characteristic cos^ (20)-dependence in the square of the bound state 
spinor wave function u{r) should be refiected in the local tunneling density of states, mea- 
sured by STM experiments in the pseudogap regime above Tc. In Fig. 2, the corresponding 
temperature evolution of |u(r)p is shown. 

Conclusions: In summary, we have studied the effect of temperature on the bound state 
wave function around impurities in dx2_y2-wave superconductors. Within the framework of 
the BdG equations, the magnitude of the binding energy is found to increase with tempera- 
ture, but stays on the order of Aq. On the other hand, the fourfold contribution to the local 
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density of states disappears gradually with increasing temperature. Therefore it appears to 
be possible to take a "snapshot" of the pseudogap by STM imaging. The present analysis 
may also be applicable to other unconventional superconductors, such as the layered organic 
superconductors k - (ET)2Cu[N(CN)2]Br and k - (ET)2Cu(NCS)2, and Sr2Ru04. 

We thank Alexander Balatsky, Ivar Martin, Bruce Normand, and Hyekyung Won for 
useful discussions. S. H. acknowledges the Zumberge foundation for financial support. 
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FIG. 1. Temperature dependence of the gap amplitude A (solid line), the bound state energy 
E (dashed line) , and the coefficient of the fourfold symmetry term in the bound state wave function 
P (dot-dashed line). The proportionality constant a has been set to a = 2/3. 
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FIG. 2. Temperature dependence of the bound state wave function, |n(r)p, localized around a 
strong-scattering impurity in a dx2_y2-wave superconductor, left: superconducting regime (T = 0), 
center: pseudogap regime (T > Tg), right: normal state (T = oo). 
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